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This paper presents the features of wind power distributions that have been analytically obtained from 
wind distribution functions. Simple equations establishing a relationship between mean power density 
and wind speed have been obtained for a given location and wind turbine (WT). Four different concepts 
relating wind power distribution functions are shown: the power transported by the wind; the 
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approximate limit to the maximum power obtained from a wind turbine, considering its parameters. 
Similarly, four different equations are obtained establishing relationships between the mean power 
density and the mean wind speed. These equations are very simple and very useful when discarding 
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Weibull distributions 
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locations for wind turbine installation. 
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1. Introduction 


WORLD wind power has increased its capacity due to the 
installation of a large number of wind turbines in recent years. As the 
rated speed of new WTs has continuously risen, it has become more 
and more important to pick the most suitable and worthwhile 
locations to install them [1]. In fact, one of the most important 
problems to be solved is whether a given location should be chosen 
or discarded for the installation of a generator system. For instance, a 
location can be valid from the point of view of wind energy potential, 
but not valid due to the lack of electrical infrastructure, whilst a not- 
so-good location according to wind speed potential can be ina much 
suitable place regarding electrical infrastructure. Due to these 
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considerations, locations close to adequate electrical lines usually 
have to be analyzed in order to investigate their energy potential, 
whilst other can be dispensable if there are signs that they will not 
cover costs of installation, taxes, etc. 

Wind power distributions have been widely investigated, 
employed and explained [2]. In many wind power studies the 
features of such distributions are used for design purposes. Both 
analytical and Monte Carlo simulation methods can be carried out, 
although they are generally used with features of wind power and 
not output power in mind. However, things can be planned from a 
different point of view, as similar distribution functions can be 
described for power, if wind distribution functions are taken into 
account, together with WT features, on the basis of data provided 
by the manufacturers. 

On the basis that the distribution function of the wind speed ina 
certain location depends just on the mean wind speed, a 
distribution function of the wind power can be obtained for a 
given WT by using its power curve. Once the wind power 
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distribution function is obtained, the mean power available is 
deduced. So as not to depend on the type of WT, this will be shown 
per unit of surface (mean power density). This process is performed 
in four different ways: 


(1) obtaining of the power from the wind [3]. 

(2) consideration of Betz’ law [4]. 

(3) consideration of realistic values [5], remembering that Betz’ 
law is an upper limit. 

(4) consideration of WT parameters such as Cut-In and Cut-Out 
wind speed, rated speed, and rated power [6]. 


The goal of this work is to obtain expressions that allow us to 
give response to questions about the mean value of the statistical 
distribution of the maximum power obtainable from the wind, 
regardless of the WT chosen, and also taking into account its 
features, when the only input value is the mean wind speed. 


2. The proposed relationships 


2.1. Relationship between wind speed and usable power from the 
wind 


The power, P’, transported by an airstream flowing with a given 
speed, U, can be calculated according to (1) as has been established 
in [7-9]. 

/ 1 3 
P= 5 Apu (1) 
where pis the air density and A the area of the airstream, measured 
in a perpendicular plane to the direction of the wind speed. 

The calculation of the mechanical power that can be extracted 
by the rotor of a horizontal axis WT (HAWT), requires Betz’ law to 
be taken into account. It states that a maximum portion of 16/27 of 
the power transported by the wind can be converted into 
mechanical power by means of such a converter. 

Thus the maximum power that can be extracted from the 
airstream is given by (2). 


8 


ye OF 
ay, 


ApU? (2) 

Eqs. (1) and (2) express the instantaneous power, P, as a 
function of the instantaneous wind speed, U. However, the wind 
speed may vary during a period of time. To consider this effect we 
are going to work with the wind speed probability distribution 
function (PDF). 

According to [10] the wind speed at a certain location may be 
represented by a Weibull distribution function with two para- 
meters called scale parameter, C, and shape parameter, k. 
Occasionally, and under certain conditions, k=2 is assumed, 
which constitutes the particular case of the Rayleigh distributions. 

Weibull PDF [11] of the wind speed at a certain location can be 
expressed as in (3). 


k U\ 1 k 
s= fele) os ee (3) 
0 U<0O 


In order to obtain the wind power PDF, a change of variables 
may be operated. As can be read in [12] and in order to obtain the 
PDF of a variable P = g(U), when g is monotonic, and where U has 
PDF fy(U), the operation given in (4) can be applied. 


Fp(P) = fulg'(P)) (4) 


t | 
g'(g-1(U)) 


where g`! denotes inverse and g’ derivative function. 


If (4) is applied to the case that is being analyzed in this paper, 
the PDF of the power may be expressed as a function of the variable 
P, i.e. P of (1) or P” of (2), considering its relationship with the wind 
speed U as in (5). 


P 


Fp(P) = fure (s 


(k/3)-1 ia 
—(P/crC?) > 
) e P 2 (0) (5) 


P<0 


where c, can be Ap/2 or 8Ap/27 depending on the case (again, P' or 
P”). 

Therefore the power transported by the wind can be 
represented by a Weibull distribution with the parameters C’ 
and k’ given in (6). 

1 k 
Ween, 3 ee 
C= zAec k 3 (6) 

And the maximum power that can be extracted by a WT can be 
represented by a Weibull distribution with the parameters C” and 
k" given in (7). 

8 k 
nw OF 3 fl D 
C= x7 ec k 3 (7) 

Moreover, as has been presented in [13-20], considering the 
features of Weibull distributions, the mean value, 4y, can be 
expressed as a function of the parameters C, k, and the Gamma 
function, such as in (8). 


1 
Umean = Huyu = cr(1 +z) (8) 


Then, the mean value of an airstream wind power PDF can be 
expressed as in (9). 


Poca er(1 | r) (9) 


And the mean value of the PDF of the maximum power that can 
be extracted by a WT can be expressed as in (10). 


1 
P” mean = Mp = cr(1 t z) (10) 


As has been explained, the shape parameter, k, of a Weibull 
distribution of wind speeds in a certain location is sometimes 
considered to be 2, giving the particular case of the Rayleigh 
distributions. Under this assumption, (8)-(10) can be expressed as 
in (11)-(13) respectively. 


Vig 


U mean =r (11) 
3 3 

P’ mean =z VTApC (12) 
ri 2 3 

P'' mean = § VTApC (13) 


and it is also possible to write the mean power as a function of the 
mean wind speed as in (14) and (15). 


3 

P' mean = =A pPUmean (14) 
16 

P” mean — 977 ^P Umean (15) 


The power can be calculated per surface unit, in order to 
establish a criterion to assess the mean power. If the value of 
p= 1.225 kg/m? for the air density (at 15 °C) is introduced in the 
previous Eqs. (16) and (17) are obtained and can be considered as 
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theoretical limits. 


Pmean _ 4173 


A mean 


(16) 


P''mean = 0.69U2 


A mean 


(17) 


The use of (16) and (17) requires using m/s as wind speed units, 
in order to obtain W/m? as power units. Eq. (16) gives the mean 
value of the PDF of the power per surface unit in a location with a 
given mean wind speed. Eq. (17) gives the maximum power that 
can be extracted, per surface unit, according to Betz’ limit. 


2.2. Relationship between wind speed and power produced by a 
HAWT 


The relationship between wind and power in a WT has been 
widely studied and depends on several items such as the type of 
WT, the air density and the temperature. 

The type of WT influences because the area swept by the rotor 
depends on the blade size. It also influences on the power 
coefficient, C,, which is the ratio between the power produced by a 
WT and the power carried by the free airstream. 

The power that a HAWT can finally extract from the wind can be 
expressed as in (18). 


P” = 5 A? pU?C (à, B) (18) 


The power coefficient is always between 0 and the Betz’ limit, 
16/27. 

In order to obtain similar results to those in the above section, 
the term Cp,,,, Will be used. This term is a constant that will 
represent a maximum value for the power coefficient. Then, the 
power P” is just a function of the third power of the wind speed 
and, therefore, it has a Weibull PDF with the parameters expressed 
in (19). 


1 k 
mW > 3 mM 
CY = SAM Page? "= (19) 


And the mean power, considering more realistic limits than 
Betz’ law, may be shown as a function of the parameter C (20). 


3 
P”' mean = 8 VITA CP imax c (20) 


It may also be expressed as a function of the mean wind speed 
(21): 


i 3 
P’ ‘mean = bP CP imax Uien (21) 


When establishing the criteria to evaluate the mean power per 
surface unit, the more realistic limit taken is shown in (22). If a 
value of 0.45 is assumed for Cp,,,, then: 


P" mean 3 
e = 0.53 U hean (22) 


2.3. Relationship between wind speed and power produced by a WT 
operating at maximum power capacity considering Cut-In and Cut- 
Out limits 


An even more realistic analysis of the wind power curves shows 
that these curves can be split in four different parts. 


e The first part corresponds to low wind speeds, and it is not 
possible, by means of any control system, to obtain power from 
the wind. The upper limit of this part is the Cut-In speed. We will 
denote this wind speed as Ucut-in. 


Table 1 
Output power in a WT depending on wind speed. 


Wind speed interval 


U < Ucut-in P=0 

Ucut-in < U < UPmax P = 3ApU?°C (à, B) 
UPmax < U < Ucut-out P=Prmax 

U <Ucut-out- P=0 


Output power 


e The second part corresponds to an interval of wind speeds, where 
maximum coefficient power can be reached (for a given variety 
of WTs). Generally, no pitch control is operated within this 
interval. 

e The third part of the curve corresponds to the maximum power 
that the WT will generate, and pitch control generally operates in 
this part. The upper limit of this part is given by a Cut-Out point 
and the corresponding wind speed will be denoted as Ucut-out- 

e The fourth part corresponds to wind speeds of values U > Ucut- 
out) and no power is generated for the sake of security. 


According to this, and in order to obtain the PDF of the 
power generated by a WT, the following considerations must be 
assumed: when U < Ucut-jn NO power can be generated. When 
Upmax < U < Ucut-out the power generated by the WT may be 
considered constant and equal to this maximum power, Pmax- 

Thus, the PDF of the wind speed may be divided into four main 
intervals which can be seen in Table 1. 

The probability of the wind speed being lower than the Cut-In 
point or greater than the Cut-Out one can be calculated as in (23) 
and (24), where probability has been denoted as Pr(-). 


Pr(U < Ucut-n) = Fu(Ucuein) = 1 — eeu /0)" (23) 


Pr(Ucut-out < U) =1- Fu(Ucut-out) = e= cur-our/O* (24) 


The probability of the power being equal to zero coincides with 
the probability of the wind speed being lower than Ucut-tn and 
higher than Ucut-out, expressed in (25). 


Pr(P =0)=1 + e~(Ucur-our/C)* _ e~ Ucut-in/C}* (25) 


The probability of the wind speed, U, being higher than Upmax 
but lower than Ucut-out is calculated in (26) according to the 
Weibull cumulative distribution function (CDF). 


Pr(UPyax, < U < Ucut-our) = Fu(Ucut-our) — Fu(UPmax) = 
— e~ (Urmax/©)" _ e-(Ucur-our/C)* (26) 


Therefore, the probability of the output power being equal to 
Pmax coincides with the probability of the wind speed being higher 
than Up,,,, but lower than Ucut-out, that has been expressed in (26). 

The probability of the output power being between the limits 
Prin and Pmax Will be taken as in Section 2.2. Pmin is the output 
power obtainable from Ucut-in, Pmin = (1/2) pACPmaxU cut-in? : 

The PDF is completed with two cases, i.e. when the output 
power is between 0 and Pin, and when it is greater than Pmax- 

Summarizing, the PDF of the output power is defined in (27). 


a+ e- (Ucut-out/C)" = e~ (Ucut-in/)") §(P) P=0 
0 O0<P< Phin 


Fo(P)=4 — (k/3cpC3)(P/c,C3) Te Pisin < P< Pmax 
(e= Urmas /O" — g~(Ucur-out/C)")§(P — Pmax) P = Prax 
0 Pmax <P 

(27) 
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Fig. 1. PDF of the output power in a WT according to (27). 


where c; = 1/2ACp,,,,,C and k are the parameters defined in (3), Pis 
the Output power, and 6 is Dirac’s delta, a generalized function, 
defined in Appendix A. 

A part of the result given in (27) has also been obtained in 
[16,19]. 

For a more detailed illustration of how the power PDF behaves, 
the result of a 20,000-sample Monte Carlo simulation is shown in 
Fig. 1. 

In order to obtain the mean output power such as in previous 
cases, it is necessary to solve (28). 


Pmean = XUP;Pr(P;) 
i=0 


where the terms P; and P, are expressed, respectively, in (29) and 
(30). 


= Pi + P2 (28) 


P, k/3 
ma (k/3)-1 _(_p 
P(k/3exC?) ( =) e (z) dP (29) 
E cC 
(e= (Urmax/ O" — @—Ucut-out/C)") P rax (30) 


In order to solve (29), the Maclaurin power series are applied 
[21], so utilizing (31), the integral shown in (29) is converted into 
(32). 


oo 


k/3 
e- Play yoo 1y aP I wat (31) 
n=0 
Pmax co 
(n+1)k 
Re | ef Ure uaa EAP (32) 
n=0 
Prin 


which can also be converted into (33). 


ke (-1)" 


Pe 
3 n(c) 


yon k/3 gp (33) 


Prin 


(n+1)k/3 


that, once solved, can be written as in (34). 


Pmax 


Pmin (34) 
(nk + k + 3) 


se k(-1)" [paea] 
Pi = 


1)k/3 
nab nl(c,C3) 


In (34) the limits can be taken into account to give (35). 


o0 k(— jr oie ES pUetEEa)/33 
Pi = ED min (35) 
n=0 ni(crC3)\ oe 


: * (nk +k-+3) 


In the case of a Rayleigh distribution (k = 2), (35) results in (36). 
pas (Phar O ~Pavin) ai 
m na an5) 


Equation (36) can be written depending on wind speed, 
resulting in (37). 
æ% 2(—1 n y5) _ ye2nts) 
P, = 5 (1) Ce ax Cut-In ) (37) 
a0 mice") (2n + 5) 


Moreover, by substituting parameter C as a function of U mean, 
(38) is now obtained. 


p, = So Mia c Ufaa — Uet) (38) 
G Mnu P (2n +5) 


Finally, replacing c; for its value, (39) may be treated as a final 
result. 


nant A 0Cp,, (Ue? = 


œx (-1) P, 
Py = max 
2 22D ust) (2n + 5) 


Ue = 


As in former sections, the mean power density is obtained in 
(40), by taking into account both terms, P4, given in (39), derived 
from (29), and P, given in (30), and finally, writing Pmax as a 
function of p, A and Cp,,,,. In this case it depends on Ucut-in, UP max: 
Ucut-out and Umean: 


Pmean _ Pi + Po 
A A 


e-UPmax/C) — e—(Ucur-out/O* 
= PCP max 


2 UPmax 
oo (=1)" mH (URH — 


prg a) (40) 


a 2DU (2n + 5) 


It can be seen that there is a dependency on parameters that 
have to be given by manufacturers in WT data sheets. In section III 
a real WT is analyzed by using these equations. 


3. Application to a real WT 


Eqs. (16), (17) and (22) show very useful relationships between 
the mean power density and the cube of the mean wind speed. 

The mean power obtained in each case has a different meaning. 
It should be remembered that the mean power that the wind 
carries is obtained in (16); the maximum mean power that may be 
obtained theoretically through a real WT is calculated in (17); and 
the maximum value of the mean power that may be obtained with 
a real WT, considering its current limits, is deduced from Eq. (22). 

A more realistic approach considers the maximum output 
power of the WT and the wind speed Cut-In and Cut-Out limits of 
the WT. In that case, and, in order to obtain an equation similar to 
those described above, the specific data of a WT have to be 
considered. 

The WT used in the simulation tests is a2 MW Vestas V80, with 
Ucut-in = 4 m/s, Ucut-out = 25 m/s. These data are given by Vestas in 
the machine data sheet. 

As the values of Ucut-in, Ucut-out and Pmax are specified, the 
results obtained are shown as a function of Umean. The values 
obtained for that WT are given in Table 2, and they are also 
displayed in Fig. 2. The results are similar to those obtained in [22]. 
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Table 2 
Values Obtained Using Eqs. (16), (17), (22) and (40). 

MWS MPD (16) MPD (17) MPD (22) MPD (40) 
4 75 44 34 30 
5 146 86 66 60 
6 253 149 114 97 
7 401 237 182 134 
8 599 353 271 168 
9 853 503 386 198 

10 1170 690 530 222 

11 1557 918 705 240 

12 2022 1192 916 253 

13 2571 1516 1164 260 

14 3211 1893 1454 262 

15 3949 2329 1789 261 

16 4792 2826 2171 257 

17 5748 3390 2604 251 

18 6823 4024 3091 243 

19 8025 4733 3635 235 

20 9360 5520 4240 226 

21 10,835 6390 4908 217 

22 12,458 7347 5643 207 

23 14,235 8395 6449 198 

24 16,174 9539 7327 189 

25 18,281 10,781 8281 180 
_ 600 
= 500 
2 
2 400 
A 
F 300 -E (16) 
o (17) 
2 200 
2 == (22) 
c 
È 109 = (40) 


t 3 & 7 Soa 45 17.39: 21. 23: 25 
Mean Wind Speed (m/s) 


Fig. 2. Mean power density as a function of the mean wind speed according to 
Eqs. (16), (17), (22) and (40). 


In Table 2, MWS means “mean wind speed”, given in m/s, and MPD 
means “mean power density”, given in W/m?. 

Although there are not so many locations with mean wind 
speeds as high as 15 m/s or more, the results are interesting in 
order to show how a WT responds to these wind speed values, 
having a maximum around 14 m/s (for that WT) and decreasing 
smoothly beyond this value. 

Table 2 shows that, for example, for a location with a mean wind 
speed of 7 m/s, there are 401 W/m? available in the wind stream; 
no more than 237W/m? may be obtained through a WT 
theoretically; no more than 182W/m? would be available 
considering real technical limits of a WT; and less than 134 W/ 
m? could be obtained through a Vestas V80-2.0, which has been 
chosen for the example. 


4. Conclusions 


The dependency between wind speed and maximum power 
obtained from a WT is described in this paper. 

In a first step an analysis of the wind power PDF has been 
carried out. It has been expressed as a function of the wind speed 
PDF, assumed to be a Weibull function. The mean power density 
has been calculated by means of an expression, which has been 
given in (16). 

The maximum power that a WT can extract from the wind was 
treated in a similar way. Betz’ law was obviously taken into 
account. The result was given in (17). 


The next step consisted of establishing a coefficient, Cp,,,,, in an 
attempt to fix a more realistic limit for the power that can be 
extracted from a WT. This limit was fixed at 0.45 for this work, but, 
technological advances will see it reaching values closer to the 
Betz’ limit in future years. The results by using this coefficient were 
given in (22). 

Furthermore, considering an even more realistic case, where 
parameters of a real WT are Ucut-in, Ucut-out, UPmax and Pmax a PDF 
was obtained for the power extracted by a WT, and the mean 
density power was obtained as a function of those parameters and 
the mean wind speed. 

As a result of this work, the following questions can be 
answered by using the result obtained: 


(1) Given the mean wind speed as a unique feature of a location, 
what is the maximum mean power that may be extracted from 
the wind? 

(2) Given the mean wind speed as a unique feature of a location, 
what is the maximum mean power that can be extracted from 
the wind, considering a given WT? 

(3) As long as the WT’s are being developed and power coefficients 
are increasing, what is the theoretical limit for the mean power 
extracted from the wind in that situation? 
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Appendix A. Appendix 
A.1. Dirac’s delta 


The Dirac delta [21] is a generalized function or density 
distribution function that can be defined as follows: 


0 t40 
TORE t=0 


faoa- faoa- 1 é>0 


Í FOS- adt = f(a) 
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